Quantum-Mechanically Induced 
Asymmetry in the Phase Diagrams of Spin-Glass Systems 
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The spin- 1/2 quantum Heisenberg spin-glass system is studied in all spatial dimensions d by 
renormalization-group theory. Strongly asymmetric phase diagrams in temperature and antifer- 
romagnetic bond probability p are obtained in dimensions d > 3. The asymmetry at high tem- 
peratures approaching the pure ferromagnetic and antiferromagnetic systems disappears as d is 
increased. However, the asymmetry at low but finite temperatures remains in all dimensions, with 
the antiferromagnetic phase receding to the ferromagnetic phase. A finite-temperature second-order 
phase boundary directly between the ferromagnetic and antiferromagnetic phases occurs in d > 6, 
resulting in a new multicritical point at its meeting with the boundaries to the paramagnetic phase. 
In d = 3, 4, 5, a paramagnetic phase reaching zero temperature intervenes asymmetrically between 
the ferromagnetic and reentrant antiferromagnetic phases. There is no spin-glass phase in any 
dimension. 
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A conspicuous finite-temperature effect of quantum 
mechanics is the critical temperature differentiation be- 
tween ferromagnetic and antiferromagnetic systems. [l], Q, 

H, lij This is a contrast to classical systems where, e.g., 
on loose-packed lattices ferromagnetic and antiferromag- 
netic systems are mapped onto each other and there- 
fore have the same critical temperature. We find that 
this quantum effect is compounded and even more ro- 
bust in spin-glass systems, which incorporate the pas- 
sage from ferromagnetism and antiferromagnetism via 
quenched disorder. 

Thus, in the present work, the phase diagrams of the 
spin- 1/2 quantum Heisenberg spin-glass systems are cal- 
culated in all dimensions d > 3. In the space of temper- 
ature T and concentration p of antiferromagnetic bonds, 
remarkably asymmetric phase diagrams are obtained, in 
very strong contrast to the corresponding classical sys- 
tems. Whereas, in the limit of d — » oo, the differ- 
entiation of the critical temperatures of the ferromag- 
netic and antiferromagnetic pure systems disappears, the 
Tp phase diagrams remain strongly asymmetric at low 
but finite temperatures, where quantum fluctuations re- 
main dominant independent of dimensionality. A di- 
rect second-order phase boundary between ferromagnetic 
and antiferromagnetic phases, also not seen in isotropic 
classical systems, is found in d > 5. In lower d, a 
paramagnetic phase intervenes between the ferromag- 
netic and antiferromagnetic systems. Our calculation is 
an approximation for hypercubic lattices and, simulta- 
neously, a lesser appro ximation for hierarchical lattices 
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The spin- 1/2 quantum Heisenberg spin-glass sys- 
tems have the Hamiltonian —(Hi. = Y^Uj) Jij s i ' s j = 
Yltij) ~f37~L{hj)i where (ij) denotes a sum over pairs of 
nearest-neighbor sites. Jij is equal to the ferromagnetic 
value of J > with probability 1 — p and to the anti- 
ferromagnetic value of —J < with probability p. We 
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the Suzuki- Takano rescal- 
Hillll to non-uniform 
systems and to length-rescaling factor b — 3, necessary 
for the a priori equivalent treatment of ferromagnetism 
and antiferromagnetism, followed by the essentially ex- 
act treatment [IjJ 1| of the quenched randomness giving 
the non-uniformity. In one dimension, 
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where the sums and products i are over every fourth spin 
along the chain, the traces are over all other spins, and 
—f3'H' is the renormalized Hamiltonian. Thus, the com- 
mutation rules are correctly accounted for within four- 
site segments, at all successive length scales in the itera- 
tions of the renormalization-group transformation. The 
trace tr is performed by quantum algebra, as given below. 

The rescaling is extended to dimensions d > 1 by bond- 
moving, namely by adding interactions resulting 
from the decimation of Eq.([T]), to obtain the renormal- 
ized interaction strength Jf, = R({Jij}), where {Jij} 
includes b d interactions of the unrenormalized system. 
The interaction constant values {Jij} are distributed 
with a quenched probability distribution V(Jij) [H, 11, 
which starts out as a double-delta function but quickly 
becomes complicated under its renormalization-group 
transformation, given by the convolution V'(J^jt) — 

J' ll/ dJ tJ V(Jij)}5(J^, - R{{Jij})). This equation ac- 
tually involves b d convolutions (for example, 729 convo- 
lutions for the d — 6 system discussed below), which 
are constituted of triplet convolutions of interactions in 
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FIG. 1: Phase diagrams of the quantum Heisenberg spin-glass systems in temperature 1/J versus antiferromagnetic bond 
concentration p for d — 3 to 10. All transitions are second-order, between the ferromagnetic (F), antiferromagnetic(AF), and 
paramagnetic (P) phases. 



series (decimation) and pairwise convolutions of interac- 
tions in parallel (bond- moving) . The quenched proba- 
bility distribution V(Jij) is kept numerically in terms of 
histograms. The number of histograms multiplicatively 
increases under rescaling, until a calculationally accept- 
able maximum is reached. After this point, the number 
of histograms is kept constant by implementing a binning 
procedure before each pairwise or triplet convolution. We 
employ a new binning procedure, in which bins are de- 
marked so as to contain equal probabilities, as opposed 
to equal interaction intervals as done previously. Starting 
from the lowest J value and moving to greater ones, his- 
tograms in each consecutive bin are combined, to interac- 
tion value J — EpjJi/Spi and imposed equal probability 
p = Y,pi = 1/nbin- In this process, histograms at the 
boundaries of bins are apportioned between the consec- 
utive bins. Thus, our calculation has 125,000 histograms 
after each decimation and 40,000 histograms after each 
pairwise bond moving. The global flows of the quenched 
probability distributions yield the phase diagrams. Anal- 
ysis of the unstable fixed points and unstable fixed distri- 
butions attracting the phase boundaries yields the order 
of the phase transitions. 

Calculations are done for the quantum Heisenberg 
spin-glass systems in integer dimensions. No finite- 
temperature phase transition occurs in d = 1,2. The 
phase diagrams for d = 3, 4, 5, 6, 8, 10 are shown in Fig[T] 
They are all strikingly asymmetric, especially in the mid- 
dle p and low-temperature (would-be spin-glass phase) 
region. In rf = 3, our calculated ratio of the critical 
temperatures of the pure antiferromagnetic and ferro- 



magnetic systems is Tq F ' /Tq = 1.48. This value is to 
be compared with the values of 1.13 found in the cubic 
lattice P, and 1.22 found in the b — 2, d — 3 hierar- 
chical lattice [1, 0| ■ This critical temperature difference 
is consistent with the lower ground-state energy of the 
antiferromagnetic system, as calculated [27[ in d = 3. 
Our calculated ratios of the antiferromagnetic and ferro- 
magnetic critical temperatures, for d = 4,5,6,8,10, de- 
crease as 1.22, 1.12, 1.07, 1.02, 1.01 respectively. On the 
other hand, it is seen that although the phase bound- 
aries leading to the pure ferromagnetic and antiferromag- 
netic critical points regain symmetry as d is increased, 
the low-temperature phase diagrams remain asymmet- 
ric. The ferromagnetic phase penetrates the antiferro- 
magnetic region at low temperatures. Thus, quantum 
fluctuations present at low temperatures favor the fer- 
romagnetic phase over the antiferromagnetic phase. In 
d > 6, a second-order phase boundary occurs directly 
between the ferromagnetic and antiferromagnetic phases, 
as is not seen in isotropic classical spin-glass systems. A 
new multicritical point occurs where all three second- 
order boundaries meet. 

The phase transition, between ordered phases, that is 
driven by quenched randomness presents a contrast to 
phase transitions between ordered phases driven by a 
systemwise uniform interaction. The latter phase transi- 
tion is obtained, at low temperatures, by driving a uni- 
form interaction that favors another ordered phase over 
the existing one. Under these conditions, essentially the 
entire system remains in one ordered phase until the 
phase transition point is reached, when essentially the 
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entire system changes over to the other phase. Through- 
out this process, the ordered domains are compact and 
have fractal dimensionality equal to spatial dimensional- 
ity, which translates to having a renormalization-group 
eigenvalue exponent of y — d, the condition for a first- 
order transition [28l |. By contrast, the phase transition 
with quenched randomness is obtained, at low tempera- 
tures, by driving quenched local interactions that favor 
the other ordered phase over the existing phase. This 
means, for example, increasing the number of random 
antiferromagnetic bonds when the system is in the ferro- 
magnetic phase. Under these conditions, the ferromag- 
netic domains avoid the random localities of antiferro- 
magnetic bonds in the system. The ferromagnetically 
ordered domains loose weight as the transition is ap- 
proached, so that the average magnetization decreases. 
At the phase transition, the ordered domains are not 
compact and have fractal dimensionality less than the 
spatial dimensionality, so that the magnetization is zero. 
This translates to the renormalization-group eigenvalue 
exponent y < d, meaning a second-order phase transi- 
tion. The converse happens when the phase transition 
is approached from the opposite side, with non-compact 
antiferromagnetic domains avoiding the random localities 
of the ferromagnetic bonds. Similarly, in another recently 
studied system with quenched randomness, second-order 
transitions between ferromagnetic and layered ordered 
phases and between antiferromagnetic and columnar or- 
dered phases, are obtained in the exact solution of clas- 
sical anisotropic spin-glasses on a hierarchical lattice. p9j 
Thus, we find that whereas phase transitions between or- 
dered phases are first order when driven by a uniform in- 
teraction, they are second order when driven by quenched 
randomness. 

In d = 3,4,5, the paramagnetic phase reaching zero 
temperature (as an extremely narrow sliver in d = 5) 
intervenes between the ferromagnetic and antiferromag- 
netic phases. In all cases, the ferromagnetic phase pen- 
etrates, reaching the high p values of 0.63 and 0.83 re- 
spectively in d = 3, where there is a zero-temperature 
paramagnetic interval, and d > 4, where there is no zero- 
temperature paramagnetic interval. The antiferromag- 
netic phase recedes at low temperatures, thereby showing 
a reentrant phase boundary [26j. 

There is no spin-glass phase, in the quantum system, in 
any dimension. The quantum version of the Sherrington- 
Kirkpatrick model [301 ] . namely the spin-1/2 quantum 
Heisenberg model with equivalent-neighbor interactions, 
with a symmetric gaussian distribution, studied from the 
high-temperature side, yields a finite-temperature phase 
transition, which has been interpreted as a transition to a 
low-temperature spin-glass phase 31] . This model should 
be similar to our studied models at p = 0.5 in the large d 
limit. Thus, we also find a finite-temperature phase tran- 
sition (Fig.l), but the low-temperature phase is explicitly 
a ferromagnetic phase with quenched bond randomness. 
The latter phase has considerable amount of short-range 
antiferromagnetic correlations, as seen in Ref.[32]. 
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TABLE I: The two-site basis states, with the corresponding 
parity (p), total spin (s), and total spin ^-component (m s ) 
quantum numbers. The state \<j>3) is obtained by spin reversal 
from \4>i). The renormalized two-site Hamiltonian —/3H(i, I) 
is diagonal in this set, with the diagonal elements of |^>i-3) 
and |04) being \j' + G' and —\J' + G' respectively. 

Appendix - The operators —f3'H'(i,l) and —f3?i(i,j) — 
(3H(j, k) — (3H(k,l) of Eq. |T]) act on two-site and four- 
site states, respectively, where at each site the spin is in 
quantum state a =j or J,. The trace tr in Eq. ([I]) is, in 
terms of matrix elements 

(u lZl \e^ H ' {l ' l) \u- Zl ) = 

Vj,W k 

(2) 

where Ui,Vj,Wk, zi,Ui,z l are single-site state variables. 
Thus, Eq. @ is the contraction of a 16 x 16 matrix into 
a 4 x 4 matrix. Basis states that are simultaneous eigen- 
states of parity (p), total spin magnitude (s), and total 
spin z-component (m s ) block-diagonalize these matrices 
and thereby make Eq. (|2|) manageable. These sets of 4 
two-site and 16 four-site eigenstates, denoted by 
and {IV'q)} respectively, are given in Tables I and II. The 
diagonal blocks are given in Tables I and III. Due to the 
microscopic randomness of the spin-glass problem, the 
four-site Hamiltonian mixes states of different parity, as 
seen in Table IV. Eq. ^ is thus rewritten as 
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TABLE II: The four-site basis states, with the corresponding 
parity (p), total spin (s), and total spin z-component (m s ) 
quantum numbers. The states |V>4,s) ; IV's), 1^13,14) are ob- 
tained by spin reversal from \ip2,i), \i>6), |^9,io), respectively. 
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TABLE III: Diagonal matrix blocks of the unrenormal- 
ized three-site Hamiltonian —/3H(i,j) — f3H(j, k) — f3H(k, I). 
The Hamiltonian being invariant under spin-reversal, the 
spin-flipped matrix elements are not shown. The addi- 
tive constant 3G at the diagonal elements is not shown. 
The interaction constants Ji,J2,J3, which are in gen- 
eral unequal due to quenched randomness, are from 
—/3Tl(i,j), —f3H(j, ft), —/3H(k, I) respectively. 
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There are only two rotation-symmetry independent el- 
ements of (^le-^'^l^p) = in Eq.©, 
which have p = p = 1,4 (thereby leading to one 
renormalized interaction constant J 1 and the addi- 
tive constant G'). From Eq. (0i||</»i) = 

(^ 1 ||^ 1 > + |(^ a ||V» 2 > + |(^3||^3> + |(^6||^> + |(^7||^> + 

i^gllV-g) - (V'gllV'lo) + KV'ioIIV'io) + 5(^16 1 1 V>16) and 
(0 4 ||04) = (MM + 2(V> 9 ||Vio> + (^ioII^io) + 

jdMIV'n) + (V»H|I^12> + 5(^1211^12) + (-01511-015), With 

(V»,||Vff> = fale-WW-PHVM-PHlWfa). From Ta- 
ble I, the renormalized interaction constant is given by 

j' = H{Wi)/(4>i\\^)). 
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